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We propose a definition of an exact lens equation without reference to a background spacetime, 
and construct the exact lens equation explicitly in the case of Schwarzschild spacetime. For the 
Schwarzschild case, we give exact expressions for the angular-diameter distance to the sources as 
well as for the magnification factor and time of arrival of the images. We compare the exact 
lens equation with the standard lens equation, derived under the thin-lens-weak-field assumption 
(where the light rays are geodesies of the background with sharp bending in the lens plane, and 
the gravitational field is weak), and verify the fact that the standard weak-field thin- lens equation 
is inadequate at small impact parameter. We show that the second-order correction to the weak- 
field thin-lens equation is inaccurate as well. Finally, we compare the exact lens equation with the 
recently proposed strong-field thin-lens equation, obtained under the assumption of straight paths 
but without the small angle approximation, i.e., with allowed large bending angles. We show that 
the strong-field thin-lens equation is remarkably accurate, even for lightrays that take several turns 
around the lens before reaching the observer. 

I. INTRODUCTION 

The phenomenon of gravitational lensing is firmly associated with the physics of a four-dimensional Lorentzian 
spacetime that satisfies the Einstein equations. Yet, it has become a common practice in the study of lensing to break 
with the basic ideas of general relativity by using the linearized Einstein equations off a fixed background, the thin 
lens approximation, and treating the bending of light as a linear phenomenon - without mention of its connection 
with the full theory. 

This point of view is very much justified by the accuracy in the comparison of contemporary observations with the 
resulting calculations, i.e., general relativity does play an essential role in lensing but the weak field approach appears 
to be quite adequate for most discussions 

However, it is now a fact that the strong field characteristics of general relativity per se are observed in nature as 
well. Black holes are possibly ubiquitous and a super-massive black hole may exist in the center of every spiral 
galaxy. Here is where the full theory of general relativity takes the leading part. In order to describe bending of light 
by black holes or in high curvature regions, it is necessary to write lens equations that respect the intrinsic nature of 
general relativity, namely: covariance and non-linearity. 

The difficulty in writing down a lens equation that respects covariance and non-linearity is very much of a conceptual 
type. In fact, even when such a lens equation is developed, it is hard to interpret. A spacetime containing a lens 
is not the superposition of two spaces, a background spacetime and a lens space. Two different spacetimes are two 
different entities, and there are an infinite number of ways of identifying them point-wise. What is the meaning of the 
angular location of a source in the absence of a lens - an idea used extensively in the thin-lens approximation? What 
are the preferred angular coordinates that give the thin-lens equation its meaning? How do we refer to the distances 
between the observer, the source and the lens in a coordinate independent manner, or what is the preferred coordinate 
distance to use? All these questions have perfectly good answers if a background spacetime is available to us and we 
are given leave to isolate the lensing action from the background. This is not so if there is no background. Without 
reference to a background, some of these questions have no answers, and some do not even make sense. Treating 
lensing phenomena strictly in the context of the full theory of relativity requires other ideas and approaches. 

We have recently introduced a proposal for a lens equation without reference to a background An exact lens 

equation on an arbitrary Lorentzian spacetime can be written down, at least in principle, since it amounts, basically, 
to finding all the light-rays that reach the eye of an observer. However, for it to be meaningful, it is necessary to 
express the equation in such a way that it can be used in an astrophysical context; it must be written or expressed in 
terms of observable quantities. To some extent, we believe that we have partially succeeded in doing that. 

As an illustration, we develop and interpret in full detail our lens equation in the case of a Schwarzschild black 
hole; explicitly working out quantities of astrophysical interest for lensing, such as the angular diameter distance and 
magnification factors. Furthermore, we use our exact lens equation to test the effectiveness of other lens equations 
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that can be written down in the case that a background is available, most notably the lens equation obtained recently 
by K. S. Virbhadra and G. F. R. Ellis [|. 

In Sec. II we discuss the idealized situation where, in principle, the null geodesic equations can be solved exactly for 
a static metric and stationary source and show, again in principle, how a set of lens equations can be constructed, 
while in Sec. Ill these ideas are then applied to the Schwarzschild black hole lensing problem. In this section, the 
important physical quantities such as the angular-diameter distance to a source and the magnification factor are 
explicitly calculated. In the subsequent sections, we compare the exact results with the thin lens calculations. 



II. THE EXACT LENS EQUATION 



We begin with a four-dimensional static spacetime (971, gah(a;°)) with local coordinates x°' and consider an observer, 
at rest in the local coordinates, on a world-line given parametrically by Xq(t), t being the observer's proper-time. 
The observer, looking out, sees null geodesies reaching him from all past null directions, 1°" . These observed directions, 
labeled by the spatial projections (orthogonal to the observer's velocity vector, v"^ = ^ccg(r)) of the null vectors, 
can be taken as the two angular coordinates of the observer's (past) celestial sphere, (01,02). The null geodesies of 
the past lightcones from the observer's worldline thus carry these labels; the points on each null geodesic are further 
labeled by the parameter along the curve, which we take to be an affine parameter s suitably normalized so that 
lav"' = 1. Thus the past lightcone of the observer has the form 

a." = X''(a:S(T),ai,a2,s) (1) 

where the functions X°(xg(T), ai , 012, s) satisfy the geodesic equation 

X^-^aX'' = 



with 

X- ^ |-X°(xg(T),ai,a2,s) 
as 

and the null condition 



gabX^X'' = 0. (2) 

The local coordinates x'^ can be chosen so that one of them, say x*^, is timelike and the remaining three a;' are 
spacelike, i = 1,2,3. In this case, the function X^ does not vanish at any point. Although we are interested in the 
past lightcone, it is more straightforward to work in terms of the future lightcone. The direction in which the lightrays 
are traced is not important in the case of interest, namely, static spacetimes. Therefore, X''^ is everywhere positive. 

This means that x^ increases monotonically with the affine parameter s. Because X^ does not vanish anywhere, 
then, by (||), at all points on a geodesic one of X* is non-zero (different i possibly in different sections of the geodesic). 
For definiteness, we label this spatial coordinate hy i — 1. From the implicit function theorem, we have that 

=Xi(xg(r),ai,a2,s) (3) 

can be inverted to obtain 

s = S{xo{T),ai,a2,x^). (4) 

The inversion will only be possible in patches, since it is possible that X^ vanishes at isolated points. This means that 
s will be, in general, a multiple- valued function of x^ . Still, this inversion allows us to reparametrize the geodesies in 
terms of just our coordinates and observation angles: 

a;° = X°{T,ai,a2TS{T,ai,a2,x^)) = X"(t, ai, 02, a;^), (5) 
x^ = x^ (6) 
= X^{t, «!, Q!2, S{t, ai,a2,x^)) = X^{t, ai,a2,x^). (7) 

with A = 2, 3. The idea is now to treat these equations as if they determine a source at the spacetime point (a;°, x^ , x^) 
in terms of the observable quantities (T,ai,a2) where we have assumed, /or the moment, that the coordinate value 
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for can be determined from observation. We will treat the source as slow moving or effectively at rest. In this case, 
Eq. (0) is defined as the lens equation. 

More specifically, we interpret this lens equation as follows. Consider a source at a spatial location a;*, emitting 
light at time x^ . We can think of the coordinate x^ as a type of radial coordinate. The remaining two coordinates 
are thus a type of angular coordinates. The emitted light arrives at the observer at a time r, in a direction (ai, 0:2). 
Eq. (0) expresses the angular location of a source at radial distance x^ in terms of the observation angles (01,0:2). 
On the other hand, Eq. (^) is an exact "time of arrival" equation; it relates the time of emission, x'^, at the radial 
location, x^, with the observer's proper time r, and arrival direction, (01,02)- 

The lens equation, Eq. (^), represents a map from the image (or observation) angles (oi, 02) to the source position 
angles, x"^. The map breaks down at locations where the determinant 

, 1 , ^ X ) /n,\ 

J(t, oi,02,x ) = det— (8) 

d(oi, 02) 

vanishes. J(r, 01,02,2;^) = defines the caustics (a three surface in four-space) of the family of past lightcones of the 
observer. The direct consequence of the break-down of the map is that multiple images of the same source can be 
observed. More specifically, often one can see an image, in direction (oi, 02), of an object that lies on a null geodesic 
before it reaches a caustic (in affine distance), while simultaneously seeing a different image, in a different direction 
(0^,02), from the same object along a different null geodesic, but, in this case, the object lies beyond the caustic in 
affine distance. The parity of the images is given by the sign of J. In background dependent calculations, J^^ is 
often interpreted as a magnification factor with respect to the "unlensed" source, but, as we have no background, this 
would not be appropriate here. (Note that J(t, oi, 02, x^) could have been calculated holding s fixed instead of fixed 
x^; the vanishing of J is independent of that choice. This follows from the general theory of Lagrangian submanifolds 
and maps.) 

The lens equation, Eq. (Q), is not yet entirely usable since it involves the (up to now) unobservable quantity x^. 
However, x^ can be expressed in terms of observable quantities through the use of the idea of distance. Though there 
are many definitions of distance in use in general relativity and astrophysics, several of them can be considered to be 
observable and we thus explore the feasibility of inferring x^ from the considerations of distance. We will investigate a 
definition of distance which is observable, namely the so-called angular-diameter distance - there being several closely 
related distance definitions 

Since we have, in principle, exact expressions for the past lightcone of the observer in terms of parameters adapted 
to the null geodesic congruence, we have a natural way of expressing the angular-diameter distance to the source in 
exact form. The angular-diameter distance is defined in terms of the infinitesimal area spanned by the observer's 
geodesic congruence at the location of the source per infinitesimal solid angle at the observer's location, namely: 



Da 



(9) 



In order to calculate the area dA^, we define two connecting vectors in the lightcone of the observer. By taking 
variations of the points on the lightcone with respect to the labels of the null geodesies in the congruence, we find the 
geodesic deviation vectors, or Jacobi fields: 

= ^, M« = (10) 

aoi ao2 

It is irrelevant to the area calculation whether s or x^ are held constant in calculating the connecting vectors - the 
difference, lying along the null tangent vectors to the geodesies, does not affect the area. 

The area dAs is the area spanned by these two vectors at the location of the source, namely, the norm of the 
wedge-product of the two vectors: 



dA, 



dai da2 



= |2 ((Ml -Ml) (M2-M2) - (Mi-M2)2)|^ rfoi do2. (11) 



If the solid angle at the observer subtended by the area dAg is given by dflo — A'(qi, 02) dai da2, where K depends 
on the choice of the coordinates, the angular-diameter distance is given by 

Dl = D\{ai,a2,T,xi)^2K-^ \{Mi-Nh) {M2-M2) - {Kh-M2f)\ . (12) 
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For sufficiently small values of x^, Eq. (^2|) is invertible, i.e., = a;"'^(ai, a2, r, D^). However, Da goes to zero at the 
caustic, so that beyond the caustic x^ is a multivalued function of Da and must be given in patches. 

The angular-diameter distance is observable, because it is related to the intrinsic luminosity L of the source and 
its apparent brightness S (total flux at the observer) via ^ 

4n{l + z)WY ^^^^ 

In principle, Eq. (|l3|), with Eq. ([l^), gives x^ implicitly as a function of observables: the angular location of the image 
(ai, 02), its redshift z(= los/ujq — 1) its apparent brightness and the intrinsic luminosity of the source. 

On the other hand, there may be situations where the intrinsic luminosity of the source is not available. In such 
cases, if there are multiple images observed, then we can make use of their relative brightness in order to estimate 
x^. For two images of the same source, lying at angles (a^^'' , oij^'') and {af\a'^^), the ratio of the fluxes S1/S2 does 
not depend on the intrinsic luminosity of the source L and can be interpreted as the relative magnification of one 
image with respect to the other one, or 

n2 r^(2) ^(2) 

= „ , (I) — m — 77- (.14^ 
Dl{a]'\4'\x^) 

Notice that {a^i\ 0^2^'') and (q;P\q;2^^) are two image directions of Eq. (|^) for a given value of the source coordinates 
{x^,x^). The inversion of Eq. (Q) is not likely to be feasible in closed form. Still, in principle, Eq. ( p^ gives x^ 
implicitly in terms of the angular location of two images and their relative brightness /ii2 = iS'i/S'2. 

(In the case of lensing at cosmological distances, it is customary to infer distances from redshifts. Even though 
we are not concerned with cosmological models in this paper, we consider redshifts as another alternative to infer x^ 
from an observable quantity. For a source on a worldline with tangent vector , emitted light of frequency Ws and 
observed frequency wq , the ratio Us / ujq is given by |^ 

LOs_ _ gab{x''{ai,a2,x^))v'^X''{ai,a2,x^) 
^0 ~ gabix^HX^ix-) 

where Eq. has been used for the source's space-time location. Equation ( |l5|) gives x^ implicitly in terms of the 
frequency of emission, the received frequency and the observed image angle. As our assumed source is at rest, its 
velocity is < = |goo(a;"(ai, as, xi))|-i/2(i, 0, 0, 0).) 

III. THE SCHWARZSCHILD CASE 
A. The Lens Equation 

We consider now the case of gravitational lensing by a Schwarzschild black hole of mass M. The line element is 

ds^ = f(r)dt^ ^— dr^ - r^(de^ + sin^^d^^) (16) 

with 

/W.l-^ (17) 

In order to take advantage of existing calculations [Q, we temporarily use coordinates (u, /) given by 

«=^(f-/ =^((-r + 2Mlog(2M-r)). (18) 

In these coordinates, the line element takes the form 
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ds^ = 2fdu^ - —dudl - -p.ide'^ +s\ii^ed4>'^). (20) 
V 2t^ 

The equations for null geodesies = in terms of an affine parameter s are equivalent to 

(21) 

(22) 
(23) 
(24) 




with the null condition x'^Xa = equivalent to 

^274 

Affi?-^CPu + e'^ + — ^ = 0. (25) 
sin 9 

The symbol (') stands for d/ds and A, B, C are three first integrals of the null geodesies, depending on the initial point 
and the initial direction. The constant C represents the freedom in the scaling of the afhne parameter s. Treating the 
observers location (uq, /qi ^Oi <^o) as the initial point, the constant B is related to the angle ip that the null geodesic 
makes with the optical axis (defined by the radial line from observer to the lens center). More precisely: 

B_ _ sinV^ 

c lo^y ^ ^ 

where is the inverse radial location of the observer. Lastly, the constant A can be related to the azimuthal angle 7 
that the direction of the lightray makes around the optical axis at the observer's location via 

^ ■ a ■ ^i^'^ fo-7\ 
— =:sm6^oSm7 , (27) 

where is the angular location of the observer. 

We can switch from the null coordinate u to the time coordinate t using i = \/2 (^ii — • Doing so and setting 

C — 1 allows Eqs. (pT[]2^ ) to be rewritten as 

sin^ ip 



sin 9q sin 7 sin ijjp 



^ = ±^^/l-( 727(1^ j'V (29) 

(30) 



loVfik)sin^9 

(9 \ sin^V (sin 00 sin 7 sin tA)^ 
I^J ^Wilo)' llf ilo)sinH ^ ' 

For those null geodesies of interest to us, i.e., those whose initial direction has a component pointing towards the 
Schwarzschild origin, the inverse radial distance / initially increases {I > 0) until the point of closest approach to the 
lens is reached, (with affine parameter value Sp). The coordinate, Z, then decreases (Z < 0) after Sp until reaching 
the source at some Sfm- However, for s < Sp and for s > Sp we have / 7^ 0, and thus the inverse radial distance I 
can be used as a parameter (in two patches, the incoming and outgoing) along the null geodesies for the purposes of 
constructing our lens equation in the manner of the previous Section. It plays the role of x^. 
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The value of / at the point of closest approach (i.e. at / = 0) is denoted Ip. If we assume, naturally, that the observer 
is located outside the last stable orbit (at r > 3M), then, for lightrays that do not cross the r = 3M orbit, the closest 
approach Ip is the smallest of the positive roots of 

(sinV')^ ll (1 - 2V2Mlp) - (1 - 2V2M Iq) = 0. (32) 

A simple analysis shows that for Iq < 3y/2M, Eq. (p2h has no positive roots unless 



sini/) > 3V2Mlo\j3{l - 2V2Mlo) (33) 

in which case there are always two positive roots, and the closest approach Ip is the smallest of them. It is 
simple to prove that Ip < (3\/2M)^^ for all ^ subject to Eq. (|3|), and that Ip — > (3\/2Af)^^ for sinip — > 



3V2Mloy3{l - 2V2Mla). See Fig. 

The Ip, which is a turning point of the coordinate I along the null geodesies, plays a major role in the following. 
First we notice that the term. 



from the u equation, can be rewritten with the role of now played by Ip, in the form 



sin^ \ ,2 J. 1 



using, from Eq. (p 



;2(1 - 2^/2Af^o) lU^ - 2V2Mlp) ' 



(35) 



We thus see that the dependence on ip is now hidden away in the inverse radial distance of closest approach Ip. This 
observation will simplify some of the calculations that follow. 

The past lightcone of an observer in coordinates = (t, I , 9, (p) in terms of the affine parameter s and initial 



directions {ip,^) could in principle be obtained from Eq. (21 - |24|) or Eq. (Bq - |3l|). This requires the integration of 



four non-linear ordinary differential equations which can not be done by quadratures. In the spirit of the Section II, 
however, we do not need the lightcone in terms of the affine parameter s, but in terms of a radial coordinate. Our 
radial coordinate is the inverse radial distance I, which is better suited for treating large distances than the standard 
r. In particular, the infinite range 3M < r < oo translates into the finite interval < Z < (3-\/2A/)~^. 

First we show how our radial coordinate I is related to the afhne length. Next, we integrate the lightcone in terms 
of I. For our purposes, it suffices to assume that by the time the lightray reaches the observer it has already passed 
by the point of closest approach Ip. 

Eq. (E3) can be integrated to obtain the affine parameter s in terms of the inverse radial distance: 



s = 2 



ll{l-2^j2Mlp) dr_ 

I, \\ Ppil - 2V2Mlp) - ^2(1 - 2V2MI') l'^ 



Pp{l - 2V2Mlp) dV 

(36) 



2V2Mlp) - ^2(1 - 2^/2M/') l'^ ' 

Equation ( ^ ) corresponds to Eq. (^ of the previous section. This represents one of the two available patches for 
s as a function of /. 

The affine parameter as a function of I is represented in Fig. || when the observer is at a distance of 30M. (We chose 
a relatively small distance in order to better appreciate the strong field effects.) The affine length goes to infinity as 
I approaches zero, in agreement with the fact that the lightray runs out to infinity. The affine length, s, is chosen 
to vanish at the observer's location. We see, in the diagram, that the affine length bulges towards the 3M radius, 
resulting in a double valued function of I. The bulge is more pronounced for lightrays that reach the observer at 
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smaller observation angles tp- The rays that come closer to the 3M radius spend more afhne time in reaching the 
observer, in agreement with the gravitational time delay. 

In the following, we obtain the past lightcone of the observer {to,lo,6o,(j>o) as a function of the inverse radial 
distance /, instead of the affine parameter s, and two angles ("0,7) specifying the direction of each null geodesic at 
the observer's location. 

Integrating Eq. ([2^ ) with Eq. ([29|), we obtain 



rl 

t = to + 2 



mi - 2V2Mlp) dl' 



k, \ /2(l - 2V2Mlp) - P{1 - 2V2MI) V2/'2(l - 2V2Mr) 

(37) 



- 2V2Mlp) dl 



\ - 2V2Mlp) - P{1 - 2V2MI) V2?'2(l - 2V2MI') 

This is the equivalent of Eq. (||) of the previous section. As a function of I, the time of arrival is double- valued (not 
so as a function of the affine parameter s); Eq. ( |37| ) represents one of the two patches. 

The integration of the angular coordinates of the lightcone is carried out in j7| . Representing the angular coordinates 
{6, (j)) in terms of the complex stereographic variables C = cot(0/2)e"^ the integration yields 

cot^ + e'^ cot QC-'"''") 
C = 1 (38) 



1 - e^T cot ^ihl^hl cot ^ 



where 0(/, /q: Ip) is 

e(Z,/o,g = ± I vr-2 



dl 



'0 ^ll{\ - 2V2Mlp) - P{1 - 2V2MI) 
dl' 



I Jll{l-2V2Mlp)-F^{l-2V2Ml')^ 



(39) 



The function Q{l,lQ,lp) depends on the observation angle "0 through Ip. The overall positive sign is taken when the 
value of the observation angle, ip, is positive, and the negative sign is taken for negative ip. This makes 0(1, Iq, Ip) an 
odd function of ip. 

Geometrically, <d{l,lo,lp) "represents" the angular position of the source relative to the optical axis, defined by 
the line between the lens and observer. The observer is considered to lie on the optical axis at 8 = tt. The relative 
angular position of a source is given by Q values between — tt and tt. Hence, Eq. (|3^) must be considered mod 27r, 
where values outside the range, —tt < & < tt, represent multiple circlings of the lens. When 8 = 0, 2tt, Att, . . . , the 
source is colinear with the lens and observer and would be observed as an Einstein Ring. For positive "ip, a value of 
8 mod 2tt between tt and represents a source located to the right of the lens, while — vr < 8 mod 27r < represents 
a source located to the left of the optical axis. 

Figure || shows a plot of 8 at fixed values of I and Iq, as a function of the image angle 0. We can see that 8 blows 
up at Ip = {iy/2M)~^ , which agrees with the fact that, as the lightrays approach the 2>M radius, they take a larger 
number of turns around the lens. Notice that 8 is a regular function of I for all / < Ip because the integrand diverges 
slowly, as {Ip — l)^^/'^. In fact, for numerical integration it turns out to be much more efficient to make a change of 
variables I — Ip — q and write 8 as 



rlp—lo 

8 = ± I TT- 2 



dq 



" ^ i2lp{\ - 2V2Mlp)q + {6V2Mlp ~ l)g2 - 2^Mq^ 
dl' 



2V2Mlp) - l'^{l ~ 2V2MI' 
In terms of the standard spherical coordinates {9,(p), Eq. (|38|) translates into: 



(40) 
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COS 9 = — cos 6*0 cos 8 + sin 6q sin Q cos 7 

sin 00 sin 00 ~ ta-n (cos (/>o sin 7 — sin 00 cos 7 cos 00 ) 

tan(p= ; — ■ ■ — — (41) 

cos 00 sin 9q + tan 8 (sin 0o sin 7 + cos 0o cos 7 cos 9q ) 

Equations (^), with Eq. ( ^0| ) are the exact lens equations for the case of a Schwarzschild spacetime and correspond 
to Eqs. (0) of the previous section. Notice that the observer is located at generic values of (0O)</'o)j which means 
that we do not choose, as is often done, the z— axis as the optical axis, the optical axis being the radial line that 
contains both the center of symmetry and the observer. This is because the spherical coordinates break down along 
the z— axis. If we chose the observer to lie along the z~ axis, then Eq. ( pl| ) reduce to cos 6* — cos 8 and tan0 = tan 7, 
and we could interpret 8 and 7 as the lens angular coordinates. However, this would result in erroneous predictions 
in the following subsections, unless we use additional care. In order to keep the remainder of this paper in the most 
transparent form, we prefer to keep the observer off the z— axis. 



B. Lensing Observables 

In this subsection, we describe the calculation of three key lensing observables from the lens equations: the angular- 
diameter distance, the relative magnifications, and the time-delay between the arrival times of two images. 

We start by exploring the angular-diameter distance, using Eqs. ( ^l|) and Eq. ( ^ to obtain an exact expression of 
the angular-diameter distance in terms of the inverse parameter I. In the next subsection, we will use the expressions 
we obtain here to explore the possibility of inferring the inverse radial distance, Z, to the source. 

First, we define the connecting vectors 

Mf.(|i.#, *#) = (0.0,*2*) (42, 



07 07 t)^ J \ 07 07 

^(didimdi\^(di mdi\ 
^ ~ \dij' d^j' di;' dtp J {dtp' ' dtp' dtp J ^ ' 

where t, 9, (p are functions of (Z, -0, 7) given by Eq. ( ^7| ) and Eq. (|4^). The partial derivatives are taken at fixed value 
of /. From the expressions above, we have 

d9 sin 8 sin 00 sin 7 ^^^^ 



97 -^/l — (cos 8 cos 00 — sin 8 sin0o 0037)^ 

90 sin 8 cos 0o + cos 8 sin 0o cos 7 98 

dtjj (cos 8 cos 00 - sin 8 sin 0o cos 7)2 dtp 

d(p sin 8 (sin 8 cos 0o + cos 8 sin 0o cos 7) 

97 1 — (cos 8 cos 00 — sin 8 sin 00 cos 7)2 

d(p sin 00 sin 7 98 

dtp 1 — (cos 8 cos 00 — sin 8 sin 0o cos 7)^ dtp 



(45) 
(46) 
(47) 

Notice that, by Eq. ( [44| ) and Eq. (|4^), the vector Mf is proportional to sin 8 for all generic values of 0o except 
for 00 = 0,7r. This means that, generically, the vector vanishes at 8 = 0, which, by Eq. (|4l]), represents source 
points (0, (p) along the optical axis. If we had chosen the optical axis as the z— axis this essential fact would not be as 
transparent. 

With the metric, Eq. (p^, we have 

Mi-Mi = -^((|^V+sin^0(^) 1 (48) 
9t V if fd9Y , ^ fdcp^ ^ 




2/2 y 97 d'^ 97 dtp 

Using Eqs. (^ - ^) the area dA^ from Eq. (^l]) can be written as 
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d^s={-j^i -TT^ - TTT + / hrr Ml • Ml #^7- (51) 



The determinant of the lens map J (Eq. (^)), 

^ _ de_d(f^ _ d(^de_ 
97 dil) 97 i9V' 

which appears in Eq. ( |5l| ) can be simphfied using Eqs. (H- 1^: 



sin 9 99 , , 

J= , =TrT- 52 

V 1 - (cos 9 cos 6*0 - sin 9 sin 6*0 cos 7)2 



The scalar product Mi • Mi can also be evaluated using Eqs. (g|-|l): 



Mi.Mi = -^ (53) 



'^As^'^(('^X-2Pf(^X\ d^d, (54) 



Thus Eq. (|51|) becomes 



The solid angle at the observer's location is dO.^ = sin^p dip dj. The angular-diameter distance is thus 

This expression for the angular-diameter distance can be simplified by showing that dt/dip can be expressed as a 
linear function of dQ /dip. We present a short, intuitive derivation of this fact here; a more formal derivation is given 
in the appendix. 

First, we notice that both connecting vectors Mf and Mf lie on the lightconc and therefore must be orthogonal to 
the null vector that is tangent to the lightrays, with components given by = (t, Z, ^, (p). The scalar product of Afj 
with is 

'dip 2/2 \^ dip ' "'^dip 



gabM^l" = ft— - {0— + sin2 ) (56) 



where 9 and cp are explicitly given by 

■ sin 9 cos ^0 + cos 9 sin 0q cos 7 sin ip 



y/1 - (cos 9 cos 00 - sin 9 sin Oq cos 7)2 (-^ _ 2 V2M/o) 
sin 0Q sin 7 sin 



(57) 



1 — (cos 9 cos 6*0 — sin 9 sin 0o cos 7)2 ^p^i _ 2y/2AIlo) 

By inserting Eq. (|5|), Eq. (|^, and Eqs. @) into Eq. (|^), we have 

, r„ /,h 1 sin?/' 99 

5abAf = — — + , ^ = ^ . (58) 



Since gabM!^^.'^ = 0, we obtain the claimed result 

9i sin V-" 99 



^ /2/2(l - 2V2M/0) 



(59) 
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Then, using Eq. ( p9| ) in Eq. ( p5| ) and Eq. (|5^), our final expressions for the area and angular-diameter distance are 



dA, 



sin 9 



de 



dip 



o , P(l - 2^/2AfO 
sm tp- 



l^il- 2V2MI0) 



dipd'y 



(60) 



and 



n2 



sin 

2/2 sin'f/' 



56 



. 2 , P{l-2V2Ml) V 
sm w ^ 



(61) 



It should be noted that the only place where Da vanishes is at sin O 0, namely, along the optical axis. The factor 



d@/dij) does not vanish anywhere and diverges at I — Ip at the same rate as the factor (1 — sin 7/) 



2 ; ;^(1-2V2M0 xl/2 



approaches zero (See Eq. (Al) in the appendix). 

From Eq. (^) and Eq. (|6l|), we see that the square of the angular-diameter distance is proportional to the Jacobian 
of the lens mapping. Because the angular diameter distance appears in the denominator of the apparent brightness, 
S", (see Eq. (p^)), a point source lying on the caustic will be infinitely magnified in the geometrical optics limit. In 
addition, our expression for the angular-diameter distance substituted into Eq. ( [l^ ) gives the relative magnifications 
for two lensed images. 

If one observes two or more images in the directions Eq- ( ^ ) can be used to define time of arrivals, \ti\. The 

subtraction of two such times defines a coordinate time delay, which can be converted into a proper time delay along 
an observer's world line. 

Among other possible candidates to useful lensing observables, which we have not concerned ourselves with, pre- 
liminary calculations suggest that the distortion of the images of small sources could be suitable for the application 
of the exact formalism as developed in this particular section. 



C. Observables and the Parameter I 



In each of our calculations of observable quantities, the non-measurable inverse radial parameter, /, plays an essential 
role. As mentioned earlier, this parameter should be eliminated in terms of observable quantities, perhaps a physical 
distance scale. 

The most direct possibility is the angular-diameter distance expression given by Eq. (|6l|). We can see that it is not 
a simple matter to invert the angular-diameter distance in order to infer I in terms of observables and lens properties. 
Nevertheless, Eq. (^) is an implicit relationship between I and the observable I?^, and can be solved numerically in 
local patches. 

A second observational way to estimate the value of /, is via Eq.(|l^. We have indicated that the ratio of the bright- 
ness of a source in two images yields an implicit equation for the source position I through the distance relationship, 
Eq. (|6l|). Hence, the parameter / may be replaced by in all calculations. 

[An alternative approach, perhaps of only academic interest, to the inverse radial distance I can be obtained 
from the redshift of the source in closed form. If we assume that the source and observer are at rest then v°, — 
|(7oo(a:"(ai,a2,xi))|-i/2(l, 0,0,0) and = |5oo(a^S|-'/'(l, 0, 0, 0), thus gahv'^X^ = |5oo|"'/'W for both and wq- 
Using the metric, we also have g^^i = l/V^ at both locations. Thus 



' 1 - 2V2MI0 \ 



Wo I 1 - 2V2MI I ' 



(62) 



which is of course the standard gravitational redshift |^ for Schwarzschild spacetime, expressed in our notation. Thus 
I can be obtained as a function of the ratio of the observed and source frequencies (and lens mass M and observers 
position Iq) by 
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D. Image and Lens Properties 



If one is interested in learning properties of an unseen lens, then the lens equation Eq. with Eq. ( |6l| ) can be used 
in conjuction with knowledge of the image properties, especially the brightness of the two main images. Because the 
brightness S of the images is proportional to an inverse power of the angular-diameter distance Da, (via Eq. (p^), 
the brighter images will be observed when the source is located near a caustic. 

We can see that the angular-diameter distance vanishes at locations where either Aff or Afj vanish. At such 
locations, neighboring rays meet. We have that M° vanishes at O — 0, which means that neighboring rays with 
the same value of ip meet along the optical axis. (In fact, all rays with that value of ip cross there, although only 
neighboring ones contribute to the intensity of the image.) This means that an observer in line with the lens and 
a source would observe an extremely bright perfect ring centered on the optical axis, at an angle ip that makes the 
function Q vanish (the Einstein ring). From Fig. ^, it can also be seen that as the image angle approaches zero, the 
source angle, 0, tends to infinity. This means that there will be an infinite number of Einstein rings, appearing in 
principle for a black hole lens, at smaller and smaller angles ~ one each time that Q passes through an integer number 
of turns, i.e., 8 = 2n7r. However, dQ/dip goes to infinity as well, which means that the additional Einstein rings get 
much dimmer as the image angle goes to zero. 

If the source is not on the optical axis, two main images will form, one on each side of the lens, because Q is an 
odd function of ■0. One image will form at large positive angle ip, whereas the opposite image forms at small negative 
angle ip, and they will have different brightness. The images at smaller angles are dimmer, because Q diverges steeply 
at the 3M radius, i.e., dQ/d4> is large. As with the Einstein rings, in addition to the two images there will be an 
infinite number of other images, dimmer and at smaller angles for a black hole lens. 

On the other hand A/f = would vanish at locations such that dQ/d^j = 0, where neighboring rays with the same 
value of 7 meet. These lie on a plane containing the source, the observer and the lens. It can be seen from Fig. ^ 
that in the case of a black hole, where the mass is contained within the 3Af radius, G is a monotonically increasing 
function of and thus there are no points where dQ/dip vanishes. Thus, we are not concerned with these caustics. 
These caustics do form in the case where a spherical lens is modeled as a uniform dust sphere with radius larger than 
3Af, and lie on the lightrays that travel through the mass, assuming the mass is transparent. 

We can now use the lens equation Eq. (^ to infer properties of an unseen dark matter or blackhole lens. Using 
Eq. ( |6l| ) with Eq. (^, and labeling the optical axis as the z— axis, we have: 

e = e(/(Af, observables),lo, ip) (64) 

where the observables might be the brightness and luminosity or (see below) the frequency ratios. For example, if 
an Einstein ring is observed at an angle Tpi, then we know both the image angle and the source's angular location, 
i.e., 9 = 0. Then Eq. (^) yields a relationship between the mass M and the inverse radial location Iq of the unseen 
blackhole. This is not particularly useful, but if another Einstein ring is observed at an angle '02, then we have two 
equations, 

= e{l{M, observables), lo,ipi) (65) 
= e(/(A//, observables), lo, ^2) (66) 

for the two unknowns Iq and M and we can infer both the location and mass of the blackhole from lensing observables. 
Notice that this exact method is necessary in order to treat multiple Einstein rings, since the standard weak- field lens 
equation yields only one Einstein ring. This fact has been observed in [ |] , and is emphasized in the following section, 
where we compare the standard lens equation with the exact approach. As for the applicability of the method, we 
refer the reader to |^ , where magnifications and angular locations of multiple Einstein rings have been calculated for 
the case of the Milky Way's galactic blackhole. 



IV. COMPARISONS OF THE EXACT AND THIN-LENS EQUATIONS 

The available approaches to lensing build on the view that the lens is a perturbation on a given background. The 
backgrounds are normally taken as Minkowski spacetime or as a Friedmann model. For definiteness, here we restrict 
to a Minkowski background. Given the flat background, one can define the locations of the source plane and the 
lens plane, at distances Ds and Dd from the observer, respectively, as shown in Fig. ^. One can further define the 
angular location of the source on the lens plane (3, and the angular location of the image on the lens plane ip. One 
can also define Dds, the distance between the lens plane and the source plane. The thin-lens approximation consists 
in considering the bending to take place at the lens plane, the lightrays being otherwise straight lines. The straight 
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lightrays bend through a bending angle a. The thin-lens approximation is justified from the point of view that the 
distances involved are much larger than the extent of the gravitational field, and is normally accompanied by the 
assumption that the image angles are small. 

In our comparisons, we consider two different approximate lens equations available. One is the standard weak-field 
lens equation 0] , and the other is a strong-field thin lens equation obtained recently in . These two lens equations 
differ essentially in the calculation of the bending angle at the lens plane. 

The standard weak-field approximation calculates the bending angle via linearized Schwarzschild. This results in 
small bending angles, which justifies a further assumption that the source angle (3 is small. Thus the standard weak- 
field thin-lens equation is also a small-angle approximation, where tan a; = smx — x for x — /3,'0,a. The weak-field 
thin-lens equation for a linearized Schwarzschild lens is 



UID. 



ds 



(67) 



In our current notation, the angular location of the source from the optical axis is 9. We need to transform (3 into a 
function of 6 in order to make a comparison with the exact lens equation. In this approximation, however, since all 
angles are small, then from Fig. 4. P/9 = Dds/Ds, thus 



We also need to express the distances in terms of the coordinates (/, ( 
assumption, we have 



(68) 



In this case, because of the small-angle 



Dds = 



V2/0' 
1 

1 



1 



72/0 \/2/ ' 

Using Eq. ( |68| ) and Eq. (71) with Eq. (|6^), we obtain the weak-field thin- lens equation in our current notation: 



(69) 
(70) 
(71) 



I 



aV2M ll 
^ [h + 1) 



This is to be compared with the lens equation, Eq. (^l]), in the case — tt. We can see that _ ^ ^ , ^ 

in the following manner. We assume that the dimensionless quantities M I = e and M Iq < e are small and make a 
Taylor series expansion of O = vr — A in terms of e: 



(72) 

that e = e reduces to Eq. (|7|) 



e{l, lo, V') « vr - A(e = 0, b*Jo, I) + e 

vyiic i^aii oiiuw uiiau uiiu z,uiuth and first orde 
further derive the next term in the expansion 



'dA' 




£2 






H 

e=0 


' 2 


9e2 



— n 



J e=0 



Ao - Ai - A2. 



(73) 



One can show that the zeroth and first order terms in the Taylor series expansion, Eq. (|73|), yield Eq. (fz^). We can 
i-To fino Tiovf foT-m in fko QVT^ona;r>n_ xhc secoud-ordcr correction, A2, is obtained to lowest order in by 



A2 = 



15 TT (M)2 



which corrects Eq. (^H) as 



1 + Iq ( 4V2M ll 157rAf2 /3 
y ^ {lo + 1) ~ 2^-^ lo + l 



I 



(74) 



(75) 



Equation ( [Tq ) translates into a second-order accurate version of the standard weak-field thin-lens equation in the 
form 



(3 = ^j- 



4:M Da 



157r M2 Dds 



DdDs^P 4 DsDji;^- 



(76) 
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The strong-field thin-lens equation [|j calculates the bending angle in exact Schwarzschild, for a lightray that 
deviates a total angle d between the incoming direction from infinity and the outgoing direction to infinity. The 
bending angle a is not necessarily small, and will increase to infinity as rays approach the 3M radius. Thus there is 
no basis for a small angle approximation of the angles (3 and a. The image angle ip is still considered small in the 
spirit of the thin lens approximation, and the lens and source spheres are still considered planes. The strong-field 
thin- lens equation is obtained in the same manner as the weak- field thin-lens equation is, from trigonometry on the 
lens diagram, but uses the exact trigonometric functions of (3 and d, rather than substitute the sines and tangents by 
their angles. The strong-field thin-lens equation is Eq. (1) in |5|, namely 

tan (3 = tan i/j — ^^-^ [tan + tan(d — ip)] (77) 



with 



2 / dl -TV (78) 

_ 2V2Mlp) - /2(1 - 2V2MI) 







Equation ( [78| ) is obtained from Eq. (10) in Q by making the identifications r and tq ■ 

Since we are interested in comparing with the exact lens equations, we need to express /3, d and Dd, Dds, Ds in 
terms of our coordinates. In the first place, since angles are not small then we have 

tan/3 = ^tan6l (79) 



and 



1 cos 9 , , 

Additionally, from Eq. (p9|), it is straightforward to see that 

d = -e(/,Zo,V')|j=io=o- (83) 
With Eqs. (^, ^ and (H), Eq. (|3) becomes 

sin 6* - cos 6* tan (6|;=;o=o + V") - 7- tan^ = (84) 

which can be manipulated to yield 

61 = V + 6|i=io=o + arcsin tan -0 cos (6|i=io=o + V")^ ■ (85) 

Equation ( p5| ) is the strong-field thin-lens equation and is to be compared with the exact lens equation d — Q. 

We can now plot our three lens equations and see how well they compare. We choose I = lo ~ ' ^^^"^^ 

corresponds to a distance of 3000Af from the lens for both the source and the observer and set M = 1 for the plots. 
We look at small image angles, which corresponds to approaches close to the 3M radius. The plots are only made for 
positive image angles. Since the lens equations are odd functions of ip, the plot for negative -0 is a reflection through 
the origin and can be inferred from Fig. |^. 

We can see in Fig. || that the weak-field equatio n (Fz^ ) deviates significantly from the exact equation both at first- 
order accuracy and at second-order accuracy, Eq. (|75|), as expected. The main reason for the second-order correction 
not to behave significantly better than the first-order accurate equation is that none of the perturbative terms blow 
up at the 3M radius, whereas the exact equation does. The perturbative approach is clearly inappropriate for the 
strong-field regime of a black hole. 
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However, the strong-field thin-lens equation Eq. ( |85[ ) agrees remarkably well with the exact equation. In our 
example, we have an error of less than one part in a thousand, for source angles as large as about six turns around 
the lens. Notice that 3000M corresponds to only 4500km for a stellar black hole, nowhere near the galactic distances 
expected for observed lensing. The thin-lens approximation appears to do very well at these relatively small distances 
and will almost certainly do much better at galactic and extragalactic distances. 



One reason why Eq. (77) performs so accurately is the tremendous distances involved. The other reason is, however, 
that the exact bending angle through infinite distances is used. Most of the effort in using Eq. (^7|) is involved in 
evaluating the integral expression of the bending angle and its derivatives, which involves the same amount of work 
required to do the exact lens equation, as shown by our comparison above, Eq. (|8^). Thus the efficiency in terms of 
computational cost is not very high. Considering that the price is to break with the covariance and non-linearity of 
general relativity, one might even find it justified to give up the thin-lens approximation in favor of the exact lens 
equation in the case of strong, spherically symmetric gravitational fields. 

A valuable, second comparison is the time delays predicted by the exact and thin lens approaches. Time delays 
are now a very important tool in astrophysics and provide a means to estimate the Hubble constant independent of 
previous methods. We will define a time delay as the observer's elapsed proper time between the arrival of a signal 
along two distinct paths connecting the source and observer. 

When the weak field, thin lens approach is applied to the Schwarzschild lens, generically, two paths connect any 
source and observer, one passing on each side of the lens. This results from the lens equation, Eq. admitting 
two values for ip for given values of /3, M, and distance parameters. The only exception is when the source, lens, 
and observer lie along the same radial line. In this case an Einstein ring is formed, and the time delay is zero by 
symmetry. The time along each of the two paths is found by integrating 

t^i ^ [ dl + (86) 



along the trajectory determined by the weak field thin lens equation. In the integration, / is the Euclidean length 
along the thin lens trajectory, and r{l) represents the Euclidean distance from the origin to a point along the path. 
The thin lens time delay, Ai*', is the difference between the two times. 

In the exact approach, there are an infinite number of pairs of geodesies connecting the source and observer because 
geodesies may wrap around the black hole many times. However, there are only two geodesies which do not circle the 
lens; these geodesies, which can be distinguished from the others, are the analogue of the thin lens paths. The exact 
time delay, Af^ , is numerically computed from Eq. ( |37| ) by taking the difference between the times elapsed along these 
two geodesies. 

Figure ^ shows the exact and weak field thin lens time delays given the same observer and source location. The 
time delays are given in geometrical units, and the /3-axis is given in radians. In this calculation, the observer is 
located at Iq = , or a distance of 3000M, and 9 = 0. Sources are located in a fixed source plane whose distance 

from the lens along the optical axis is also SOOOAf, as in the standard thin lens picture. The time delays are plotted 
against /3, the source angle used in the thin lens equation. 

At /3 = 0, Einstein rings are formed, and each time delay will be zero. As /3 increases, the thin lens time delay slightly 
overestimates the true value, as seen in the figure. This represents the general behavior for a single Schwarzschild 
lens at larger distances, but the effect remains fairly small. For a single Schwarzschild lens with a mass of 2.5 x 10^^ 
solar masses, at redshift z = 0.5 with a source at z = 1.0 and f] = 0.5", the overall exact time delay is close to 400 
days, while the error introduced by using the thin lens approximation is less than an hour. 

The time delays computed in the thin lens approximation remain quite accurate in our comparisons because the 
light rays do not feel a strong gravitational field along their trajectories. This is why the thin lens fails so dramatically 
in the observation angle calculation presented in Fig. ^, while fairing well in the time delay comparisons. 

Despite these results, it may be incorrect to assume that the thin lens method of computing time delays remains 
accurate for more complex lensing situations. The first of two possible problems is that the weak field assumption 
may not always be valid along an observed null geodesic. It may be very difficult to identify lens candidates whose 
null geodesies have undergone interactions with strong gravitational fields, however, if such candidates are found, the 
thin lens methodology will most likely prove quite inaccurate. Secondly, it is not known how accurate the time delay 
computations will be in more complex lensing scenarios, especially when the geodesies are bent in multiple lens planes. 

V. DISCUSSION 

To date, virtually all applications of gravitational lensing have utilized the thin lens approximation. Although the 
thin lens method has proven a quick and useful tool, it can not be applied to high curvature regions. The failure 
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of the next order correction to the thin lens equation, Eq. (|7^), as contrasted with the abihty of the strong field 
thin lens equation 1^ , to capture the divergence of the observation angle Q in Fig. |^ emphasizes the point that some 
combination of exact methods are required in such situations. 

We are putting forward the idea that exact gravitational lensing may not be just a purely academic exercise. In 
fact, we see that virtually all of the observationally relevant quantities can be determined analytically in the exact 
method with relative ease for a Schwarzschild spacetime. In this paper, wc found straight-forward expressions for the 
time delays, observation angles, angular-diameter distances, and relative magnifications. These analytic expressions 
can be used in comparison calculations or model building with great computational prowess, power, or time. 

The limited testing of the thin-lens methodology in this paper has indicated that there are regions (possibly not 
yet observable) where the thin-lens method fails in a highly symmetric case. In ongoing calculations, we are exploring 
the accuracy of the thin lens approximation under a broad range of conditions, including those in which there may 
be stronger fields or multiple lenses present. These calculations will require a combination of analytic and numerical 
results. A precedent for this kind of work was set by Ranch and Blandford ||l^, who studied the null geodesic equations 
(or the exact lens equations in our terminology) for Kerr spacetime and found the caustics of the lightcone. Our sense 
from these comparisons is that the error in the thin lens method for the time delays and observation angles will, in 
some cases, be appreciable [^. 

We are also interested in an issue regarding the comparison of the sizes of two different types of corrections to the 
thin lens equations. On one hand, within the framework of the thin lens methodology, there are corrections due to the 
structure of the mass distribution (of the lens) over that of the monopole moment. On the other hand, even for the 
monopole case, there are differences between the predictions of the exact lens equations and the thin lens equations. 
The issue is whether the sizes of these corrections are comparable. 
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An alternative derivation of Eq. (p9|), by direct computation, is obtained by taking the derivative of O and of t with 

respect to ip. Because both t and 8 depend on ijj only through the point of closest approach Ip, we actually need to 
calculate dt/dlp and dQ/dlp. We have 
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APPENDIX A: 






(Al) 



On the other hand we also have 



dL 





Ppil - 2V2Mlp) ^ 

Ipjl- 3V2Mlp) 



dl' 



^l^{l-2V2Mlp)(^ll{l - 2V2Mlp) - ^'2(1 - 2V2Af?')) 



3/2 



+ 



y2Z'2(l - 2V2MI')) \ l^il ~ 2V2Mlp) - - 2^/2Ml') 





(A2) 
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Clearly the last term inside the braces under the limit sign in (^3) can be substituted with 



2/2(1 - 2V2Mlp)) Jllil - 2V2Mlp) - - 2^Ml') 



(A3) 



l'=l„—e 



since the difference vanishes in the limit e ^ 0. We can thus see that dt/dlp and dQ/dlp are proportional to each 
other via 



dL 



1 



99 



P \/2/2(l - 2^/2A/^p) 



(A4) 



This implies Eq.(& 
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FIG. 1. A plot of the function F ■- 1 - sin^ f / (ll f (h)) as a function of / for a sequence of values of sin-i/) including 
smij) = 1 . Notice that for large enough ^ there are always two positive roots, while for small enough i/) there are no positive 
roots. There is a critical value of for which there is only one positive root. The allowed values of I are such that F is positive, 
and such that Z = is included. Thus, allowed values of I range from to the smallest positive root. The smallest positive root 
is denoted Ip and represents the point of closest approach to the lens. Notice that Ip increases with decreasing xl), reaching the 
critical value (3\/2M)~^ at the critical value of iji. 
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FIG. 2. The afiine length s as a function of the inverse radial distance 1 for three null geodesies that reach the observer from 
positive directions We have set the observer at lo = 0.1/3\/2M, the mass is M = 1, and the constant C has the value 1 
for the three null geodesies. The three null geodesies are labeled according to the value of sin^p. Smaller positive angles ip 
reach closer to the 3M radius at their respective point of closest approach Ip. For this observer the smallest image angle is at 
sintp = 0.16733201, at which the closest approach reaches the 3M radius. We see that the affine length is double valued as a 
function of the inverse radial distance. 
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-2.0 -1.0 0.0 1.0 2.0 

FIG. 3. The function Q{1, lojtp) at fixed values of / and lo. The angles on both a^xis are measured in radians. We have chosen 
I = lo = 0.1/3-\/2Af with M = 1 for this figure. This plot represents the exact lens equation in the case that the a— axis is 
chosen as the optical axis, defined by the observer and the lens. 
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FIG. 4. The lens diagram for the case in which a background is available. The quantities P,ip,a,Ds,Dd,Dds represent the 
Euclidean angles and angular diameter distances in the background space-time shown. The quantity 9 represents the angular 
position of the source with respect to the z — axis in Schwarzschild coordinates if Schwarzschild spacetime is thought of as 

superimposed on a flat background. 
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FIG. 5. A comparison of the available lens equations in the literature. The vertical axis represents the angular position of 
the source in radians. The horizontal axis represents the observed angle ip in radians. We have taken the observer and the 
source at the same distance from the lens, r = 3000M. Here TL:=thin lens. The standard thin-lens-weak-field approximation is 

clearly inaccurate, whereas the second-order correction to it docs not do significantly better. However, the thin-lons-strong-field 
approximation appears to be remarkably good, since at this resolution it agrees with the exact lens equation up to one part in 
a thousand. 
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